Planck Scale Physics and Bogoliubov Spaces in a Bose— Einstein Condensate 
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We analyze the consequences caused by a deformed relation, suggested in several quantum gravity 
models, upon a bosonic gas. Concerning the ground state of the Bogoliubov space of this system, 
we deduce the corrections in the pressure, the speed of sound, and the corresponding healing length. 
Indeed, we prove that the corrections in the relevant thermodynamic properties associated to the 
ground state, defines a non trivial function of the density of particles and the deformation parame- 
ters, allowing us to constrain, in principle, the form of the energy-momentum dispersion relation. In 
addition, we calculate the condensation temperature associated to the non-interacting system, and 
show that this fact could be used also, to infer representative bounds for the deformation parameters, 
under typical laboratory conditions. 
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I. INTRODUCTION 

The possibility of a deformation in the dispersion rela- 
tion of microscopic particles, appears in connection with 
the quest for a quantum theory of gravity [H-Q. In 
some schemes, the possibility that the space-time could 
be quantized, can be characterized, from a phenomeno- 
logical point of view, as a modification in the dispersion 
relation of microscopic particles [1, H, Q (and references 
therein). A modified dispersion relation emerges as an 
adequate tool in the search for phenomenological conse- 
quences caused by this type of quantum gravity mod- 
els. Nevertheless, the most difficult aspect in the search 
of experimental hints relevant for the quantum-gravity 
problem is the smallness of the involved effects [1, 0]. 
If this kind of deformations are characterized by some 
Planck scale, then the quantum gravity effects become 
very small 0, 1. In the non-relativistic limit, the dis- 
persion relation can be expressed as follows @ 



E : 
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in units where the speed of light c = 1, being Mp 
(~ 1.2 X Ky^^eV) the Planck mass. The three parame- 
ters ^1,^2, and ^3, are model dependent and should 
take positive or negative values close to 1. There are some 
evidence within the formalism of Loop quantum gravity 
[^-Q that indicates a non-zero values for the three pa- 
rameters, a, a, a, and particulary 0,0 that produces a 
linear-momentum term in the non-relativistic limit. Un- 
fortunately, as is usual in quantum gravity phenomenol- 
ogy, the possible bounds associated with the deformation 
parameters, open a wide range of possible magnitudes, 
which is translated to a significant challenge. 
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In a previous report [37[, we were able to prove, that 
the condensation temperature associated to the non- 
interacting bosonic gas, trapped in a generic three dimen- 
sional power law potential, is corrected as a consequence 
of the deformation in the dispersion relation. Moreover, 
this correction described as a non-trivial function of the 
number of particles and the shape associated to the corre- 
sponding trap it could provide representative bounds for 
the deformation parameter £^1. We have proved that the 
deformation parameter ^ can be bounded, under typical 
conditions, from < 10^ to j^j < 10^, by using differ- 
ent classes of trapping potentials in the thermodynamic 
limit. In the case of a harmonic oscillator-type potential, 
we have obtained a bound up to 1^1 < 10'*. 

In references [1, Q it was suggested the use of ultra- 
precise cold-atom-recoil experiments to constrain the 
form of the energy-momentum dispersion relation in the 
non-relativistic limit. There, the bound associated to 

is at least, four orders of magnitude smaller than the 
bound associated to a Bose-Einstein condensate trapped 
in a harmonic oscillator obtained in reference [s^l • How- 
ever, the results obtained in reference [13], suggest that 
the many body contributions would allow to improve, 
in principle, the bounds associated to £,1. These facts 
open the possibility to explore the corrections caused by 
an anomalous dispersion relation, on the relevant ther- 
modynamic properties associated to N-body systems, for 
instance, within the Bogolioubov formalism associated to 
Bose-Einstein condensates. 

The phenomenon of Bose-Einstein condensation, from 
the theoretical and experimental point of view, has pro- 
duced an enormous amount of interesting publications 
associated to this topic [l0l - l30| (and references therein). 
Among the issues addressed we may find its possible 
use as tools in the search of quantum-gravity manifesta- 
tions, for instance, in the context of Lorentz violation or 
to provide phe nomenological constrains on Planck-scale 
physics [UHll. 

For this purpose, we define the next modified N-body 
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Hamiltonian 

H = 



n 

k=0 

2V 



2m 



atata^, -a^ r, 
k=a p=o 9=0 



(2) 



where Un 



being a the s-wave scattering length, 



mc 
2M„ ' 



mc is the rest mass, and the term ahk with a = 
in ordinary units, is the leading order modification in 
expression ([T]). The system under study will be a Bose- 
Einstein gas enclosed in a container of volume V, having 
a deformed dispersion relation as a single particle energy 
spectrum. Additionally, the interaction between two par- 
ticles will be assumed to be dominated by s-scattering, 
i.e., the temperature of the system is very low {ka << 1, 
where k and a are the wave vector and the scattering 
length, respectively) (43| . 

The main goal of this work is to analyze the corrections 
caused by a deformed dispersion relation upon the prop- 
erties associated to the ground state of a Bose-Einstein 
condensate, within the Bogolioubov formalism. In other 
words, we will show that the many-body contributions 
associated to this system, could be used, in principle, 
to establish some meaningful bounds for the deforma- 
tion parameter ^i. Finally, we calculate the shift on the 
condensation temperature for non-interacting systems, 
caused by a deformed dispersion relation, and show that 
this alternative procedure can be used also to provide 
representative bounds for ^i. 



This Hamiltonian can be diagonalized introducing the 
Bogoliubov transformations (43| 
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hi 



4 + 7fca_fe 
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These two operators fulfill the same algebra related to a^: 
and at, i.e., they are also bosonic operators. In this last 
expression the following definitions have been introduced 
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The final form for our Hamiltonian is 
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The last summation diverges, a result already known [45l . 
l46j . and this divergence disappears introducing the so- 
called pseudo-potential method, which implies that we 
must perform the following substitution [46| 



II. MODIFIED BOGOLIUBOV SPACES 

In this section let us calculate the corrections on 
some relevant thermodynamical properties, caused by 
the deformation parameter a in a Bose-Einstein con- 
densate within the Bogolioubov formalism (see for in- 
stance, [H, H^). After some assumptions, the second 
right hand term in equation ([2]), very close to the tem- 
perature T = 0, can be re-expressed as follows [3, EH] 



fe=0 P=0 q=0 
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With these approximations the iV-body Hamiltonian has 
the following structure 
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Finally, this last summation will be approximated by an 
integral. It is noteworthy to mention that the usual in- 
tegral, has as a lower limit the value k = [l^, l46j . 
However, the modified integral given here, does not have 
as a lower limit the value fc = 0. In other words. 
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where 

Additionally 
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note that in the case a = 0, we recover the usual result 
[46j , but corrected by the contributions of the rest energy. 
With these conditions we deduce the final structure of the 
A^-body Hamiltonian 



H = Eo 



(15) 



In this last expression denotes the energy of the 
ground state of the corresponding Bogoliubov space (43| , 
corrected by the contributions of the deformation param- 
eter a and the rest energy 
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a^T^ 0.05 



(0.05 + ln5) 



On the other hand, we have that the modified-energy of 
the Bogoliubov excitations (Ek) is given by 



Ek 



2UaN 



(17) 



Rescaling the energy by the rest mass and in the long- 
wavelength limit fc ^ 0, expression p7)) becomes 



Ek^ 
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167ra7V 
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(18) 



On the other hand, at the high-energy limit, fc — >■ cxi, 
expression (fT7| becomes 



Ek„. — 



2m 



thk 



mV ' 



(19) 



which basically corresponds to the deformed dispersion 
relation, expression ([T]), corrected by the mean field con- 
tributions. If we set a = in equations and ([TO]), 
we recover the usual result [i^. From equations 
and p9)) . we are able to define the corresponding healing 
length 



TTn.n, V 



1 — a- 



h\/ Snan 



(20) 



if we set a = 0, we recover the usual result [42| - |44l . |46| . 
Under typical conditions the healing length is approxi- 
mately (or bigger than) one micrometer Al]. The order 
of magnitude associated to the correction caused by the 
deformation term is of order 10 ^'^m! for |^i| < 1 in the 
case of igK for typical densities 10^^ — 10^^ atoms per 
cm^. Additionally, these facts allow us to bound the 
deformation parameter up to |^i| < 2.4 x 10~^ for a 
typical healing length of order 10^^ m. Notice that more 
diluted systems could be used to improve the bound 



associated to ^i, and also obtain a possible measure of 
the deformation parameter ^i. 



and 

TOi^^) associated to the 
ground state of the Bogoliubov space become, respec- 
tively 



On the other hand, the pressure (Pq 
speed of sound {vs — 
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The detection of the speed of sound in condensates has 
already a long history [l^ . Notice that the possibil- 
ity of detecting the term depending upon the deformation 
parameter a (5vf) requires that, if A{vs) is the experi- 
mental error, then Avs < \Svf\. In our case this entails 



A(vs) < 



a- 



Sa^hN 
irmV 



V 20 



(23) 



Under these conditions an experimental uncertainty 
A(vs) of order 10~^-y/|^i| could be tuned, in principle, 
below Planck-scale induced speed of sound, under typi- 
cal conditions. In other words, high precision measure- 
ments are required. Here, it is noteworthy to mention 
that the accuracy in the speed of sound measurements is 
not reported in the literature, at least in the literature 
known by the author. The speed of sound in a conden- 
sate is typically of order 10~^ms~^ [2^, this fact allows 
us to infer a bound for the deformation parameter up to 
161 < 6 X 10^, for tymcal densities of order 10^^ - lO^'^ 
atoms per cm^ [Tol. l47j. Conversely, the case < 1, 
lead to densities of order 10"'^* atoms per cm^, assuming 
a speed of sound of order 10~'^ms~^, that is, three or- 
ders of magnitude bigger than typical densities. Notice 
also that, the correction upon the speed of sound caused 
by a, is a non-trivial function of the density n, which 
is proportional to the product anlnn^^/^, this relation 
between the deformation parameter a and the density 
n, suggests that denser systems would allow improve, in 
principle, the bound associated to 6- 



III. CONDENSATION TEMPERATURE FOR 
THE CASE Uo = 

Let us analyze the condensation temperature associ- 
ated to an uniform Bose-Einstein gas in a box, in order 
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to obtain representative bounds associated with the de- 
formation parameter ^i, within the semiclassical approx- 
imation [lO, Hi- Let us take into account the leading 
order modification in expression ([TJ, or equivalently, if 
we set a = in expression (flQ)) . this leads to the follow- 
ing dispersion relation between energy and momentum 



Er, — 



2m 



ap. 



(24) 



The single-particle phase-space distribution for bosons, 
can be written as Qii Sjl 



n{r,p) 



1 



(25) 



where /3 = 1/kT, k is the Boltzmann constant, and T the 
temperature. Additionally, fi is the chemical potential. 
In the semiclassical approximation, the number of parti- 
cles in the 3-dimensional space obeys the normalization 
condition [1^ El] 



N = 



1 



(27rft)3 



(fif(fipn{r,p), 



(26) 



Thus, after some algebraic manipulation, we are able to 
express the total number of particles associated with our 
system as follows 



N = No + V 



93/2{z) 



m?KT 



(27) 



where V is the volume of the container, Nq is the num- 
ber of particles in the ground state, and z = e^^ is the 
fugacity 42]. The function qyjz) is the so-called Bose- 
Einstein function defined by [42[ 



1 



Notice that the Bose-Einstein function diverges for z = 1 
when v < I The behavior of the Bose-Einstein 

functions and its relation with the value of the chemical 
potential is usually used as a criterion of condensation 
[27I [4^ . Apparently, Bose-Einstein condensation in the 
thermodynamic limit is not possible for this system, due 
to the divergent behavior of the Bose-Einstein function 
gi{z), when /x = 0, as can be seen from expression pT]). 
In the usual case, a = 0, we have only the term propor- 
tional to <73/2(-z) in expression (j27p . in this situation, the 
condensation is possible in three dimensions for a bosonic 
gas in a box, and if we assume that the ground state en- 
ergy for this system is zero then, we are able to set the 
value of /i = safely without problems of divergence, 
and 33/2(1) is just C(3/2) at the condensation tempera- 
ture. Nevertheless, the minimum value associated to the 
energy for a modified bosonic gas in a box is not zero, 
in fact this value is ma^/2. These assumptions allow us 
to express the number of particles at the condensation 



temperature as follows (assuming A^o ~ at the conden- 
sation temperature) 



( 



N _ 



53/2 ( 



(29) 



By using the properties of the Bose-Einstein functions 
when Pcuna^ /2 — >■ [l^], we may now express equation 
as follows 



V ~ \ 2^ 



3/2 iti^kT 



V ma^ 



(30) 



where Tc is the condensation temperature. If we set 
a = 0, we recover the usual result for the condensation 
temperature associated to a bosonic gas trapped in a box 
in the thermodynamic limit [42i, Eij 



K.Tn = 



2T:h^ ( N 



■Wc(3/2) 



2/3 



(31) 



At this point it is noteworthy to mention that the ther- 
modynamic limit is defined as — >■ 00, V 00, keeping 
N/V constant [l^. El]. With this definition, the con- 
densation temperature is well defined. From expressions 
([501) and (PT|) . we obtain the shift on the condensation 
temperature caused by the deformation parameter a in 
function of the number of particles 



T, - To AT, 



2™ 



Tn 



3h\ N J 



-1/3 



(32) 



1-ln 



N \2/3 



(ma)2 VyC(3) 



(28) Clearly, in the limit a 0; ^ ^ 0. The current 
high precision experiments for igK, show that the shift in 



the condensation temperature respect to the ideal result, 
caused by the interactions among the constituents of the 
gas is about 5 x 10~^ with a 1% of error [s^. These facts, 
allow to compare the results given here, in order to infer 
representative bounds for the deformation parameter ^1. 
Under typical conditions, the particle density n = N/V, 
varies from 10^'^ — 10^^ atoms per cwP' [13, SS]- Using 
these data, we are able to bound the deformation param- 
eter up to 1^1 1 < 5.2, for n « 10^'^ atoms per cm^ , that is, 
two orders of magnitude less than in 37[ , and of the same 
order of magnitude that in which is noticeable. In 
addition, the correction in the condensation temperature 
caused by a, which is proportional to n~^/'^(l — Inn^/'^), 
can be enlarged in systems with sufficiently small n, but 
where the thermodynamic limit is still valid. These facts 
suggest that more diluted systems could be used, in prin- 
ciple, to improve the bounds associated to ^1 or even to 
obtain a direct measure of such parameter. 
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IV. CONCLUSIONS 

We have analyzed a modified bosonic gas in a box 
within the Bogolioubov formalism. We have calculated 
the corrections on the ground state energy of this system, 
and consequently the corrections in the pressure and the 
associated speed of sound. Additionally, we have ana- 
lyzed the shift in the condensation temperature caused 
by a deformed dispersion relation. Concerning the anal- 
ysis of the associated speed of sound, this allows us to 
bound the deformation parameter up to |^i| < 10"^ as- 
suming a speed of sound of order 10^'^ms^^. Neverthe- 
less, expression ([^5)) suggests that denser systems would 
allow to improve the bounds on , but where the dilute- 
ness condition, n|a|^ << 1, must be fulfilled. By using 
equation (P^ . we obtain for different values of the den- 
sity, n ~ 10^*^ atoms per m"^ a bound up to |^i| < 10^, 
for n ~ 10^'^ atoms per m^; |^i| < 0.7, and for n ~ 10^^ 
atoms per m'^ ; |^i| < 10~^, assuming a typical speed of 
sound of order lO^^ms^^. Additionally, the corrections 
on the healing length lead to a bound up to |^i| < 10"^ 
assuming a healing length of order lO^^m. 

On the other hand, the shift on the condensation tem- 
perature caused by the deformation parameter, leads to 
a bound up to |^i| < 5.7, under typical laboratory con- 
ditions. This bound could be improved, even more in 
diluted systems, that is, the opposite behavior, if one 
wants to infer a bound on by using the speed of sound 
or the healing length corrections. 



Clearly, the results given in this work, must be ex- 
tended to a more realistic situation. First, from the ex- 
perimental point of view, there is no condensates in a 
box. Usually, the confinement of the condensate can be 
obtained by using harmonic traps, among others H] 
and the use of generic potentials, under certain condi- 
tions, could be used to constrain Planck scale physics in 
this context, by changing the shape of the trap [3^ W^- 
In other words, the analysis of a system trapped in 
generic potential within the Bogolioubov formalism de- 
serves further investigation, and will be presented else- 
where [4(|. Finally, we must add that the measure of 
the corrections caused by the deformation parameter ^i, 
could be out of the current technology. Nevertheless, 
it is remarkable that the many-body contributions in a 
Bose-Einstein condensate, open the possibility of plan- 
ning specific scenarios that could be used, in principle, 
to obtain a possible measure of the effects caused by the 
quantum structure of space-time. 
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